Loss-induced phase separation and pairing for 3-species atomic lattice fermions 



A. Privitera, 1 ' 2 ' 3 I. Titvinidze, 1 S.-Y. Chang, 4 ' 5 S. Diehl, 4 A. J. Daley, 4 and W. Hofstetter 1 

1 Institut fur Theoretische Physik, Johann Wolfgang Goethe- Universit at, 60438 Frankfurt am Main, Germany 
2 Dipartimento di Fisica, Universita di Roma La Sapienza, 1-00185 Roma, Italy 
3 Democritos National Simulation Center, Consiglio Nazionale delle Ricerche, 
Istituto Officina dei Materiali (CNR-IOM) and International School for Advanced Studies (SISSA), 1-34136 Trieste, Italy 
* Institute for Quantum Optics and Quantum Information of the Austrian Academy of Sciences, 
and Institute for Theoretical Physics, University of Innsbruck, A-6020 Innsbruck, Austria 
5 Department of Physics, The Ohio State University, Columbus, OH 43210, USA 

(Dated: December 22, 2011) 

We study the physics of a three-component Fermi gas in an optical lattice, in the presence of a 
strong three-body constraint arising due to three-body loss. Using analytical and numerical tech- 
niques, we show that an atomic color superfluid phase is formed in this system and undergoes phase 
separation between unpaired fermions and superfluid pairs. This phase separation survives well 
above the critical temperature, giving a clear experimental signature of the three-body constraint. 



PACS numbers: 03.75.Lm, 37.10.Jk, 67.85.Lm 

The recent realization of multi-component degener- 
ate Fermi mixtures in experiments [TJ |2] offers exciting 
new opportunities to observe physics beyond that con- 
ventionally studied for two-component fermions in the 
solid state. Three-species mixtures can, e.g., give rise 
to a competition between bound three-body trions and 
an atomic color superfluid (ACS) phase, where unpaired 
fermions coexist with superfluid (SF) pairs of different 
components [3j [4] . The current limitation for observing 
this physics in three-component Lithium gases is the high 
rate of three-body loss due to a three-body Efimov reso- 
nance in this system [U[5]. However, loading the gas into 
an optical lattice could be used to suppress losses, as a 
large rate of onsite three-body loss can prevent coherent 
tunneling processes from populating any site with three 
particles [6j [7] , as was demonstrated for the case of two- 
body losses in Feshbach molecules [8 j. For three-body 
loss, this mechanism would provide an effective three- 
body hard-core constraint [10] , suppressing the ac- 
tual loss events. 

Here we address the key question of how this loss- 
induced constraint affects the many-body physics in a 
lattice in 2D and 3D. We focus on a regime of medium- 
to-strong attractive interactions with pronounced asym- 
metries in the interactions between different species, cor- 
responding to the regime of strong losses in current ex- 
periments with 6 Li [5 . On a lattice, we find that the loss- 
induced constraint not only inhibits trion formation, but 
also drives the homogeneous ACS phase towards phase 
separation of SF pairs and unpaired fermions (see Fig. [I]). 
In contrast to two-species Fermi mixtures with attractive 
interactions [TT] . phase separation takes place here even 
for a number-balanced mixture. Moreover, this separa- 
tion persists well above the ACS critical temperature, 
leaving a strong many-body signature of the constraint 
accessible at current experimental temperatures. 

Below we discuss the model and establish a simple ana- 
lytical picture for the phase separation based on deriving 
an effective low-energy Hamiltonian at strong-coupling. 




FIG. 1: (Color online) Sketch of the phases (a) without and 
(b) with the three-body constraint induced by losses. The tri- 
onic phase expected in the strong-coupling regime is replaced 
by a phase-separated ACS. 



We then analyse the system quantitatively for realistic 
parameters using Dynamical Mean-field Theory (DMFT) 
and Variational Monte Carlo (VMC) techniques, and 
present experimental signatures. 

Model - A three-component Fermi gas loaded into the 
lowest Bloch band of an optical lattice can be described 
by the Hamiltonian (h = 1) 

(i,j),cr cr,i a<a' ,i 

(i) 

where c\ a is a fermionic creation operator for an atom 
on site i in internal state a = 1,2,3, J is the tun- 
nelling amplitude between nearest neighbour sites (z, j), 
U(j jCr ' are the pairwise onsite interaction energies, and 
h a ^ = c^i&cr^. The chemical potentials \i a fix the filling 
factors n a . This model assumes that J, <C where 
v is the energy separation to higher Bloch bands. In the 
presence of strong on-site three-body loss with a rate 73, 
tunnelling processes populating sites with three atoms 
are suppressed as J 2 /73 [6 • We thus obtain a model 
with an effective three-body constraint, n^in^2^i,3 = 0, 
as discussed in Refs. [6j|9j[T0]. 

Parameters for Lithium - The requirement 73 ^> J 
is fulfilled if we take 6 Li atoms at a magnetic field of 
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500 Gauss as realized in Refs. [1 j. From experimentally 
measured scattering lengths and loss rate constants PQ, 
we obtain for an isotropic cubic lattice in 3D of depth 
SEr in each dimension (where Er is the recoil energy), 
that U la ~ -41.6 J, J7i >3 « -25.1JC/ 2 ,3 ~ -6.2J, and 
73 ~ 100 J [TO] . Variations in the lattice depth corre- 
spond to rescaling the values of the interactions given 
above by a common factor A, i.e. — >• \U a ^'. 

Throughout this work we focus on a globally balanced 
mixture, n a = n, in the medium-to-strong coupling 
regime (A G [0.3,1]). For comparison of 3D and 2D re- 
sults, we keep U a ^ a f /W constant, where W = AD J is the 
width of the lowest Bloch band in D dimensions. 

Simple picture of pairing and phase separation - To es- 
tablish an analytical understanding, we first choose A = 1 
and n a = n = 1/6. In this regime the system exhibits 
clear energy scale separation, (i) being interaction domi- 
nated \U aj(T f\/J ^> 1 and (ii) having strongly anisotropic 
interactions. 

We can see how this scale separation leads to pairing 
and phase separation by generalizing a recently devel- 
oped technique for treating onsite constrained models to 
fermionic systems [9]. The constrained fermionic theory 
is exactly mapped onto a theory in single fermion oper- 
ators t ay i and a vector boson b a ^ made up of fermions 
of species a + l,cr + 2 on site i (with cyclic property 
<j + 3 = a) [12] without additional constraints. Due to 
the balanced mixture and strong scale separation, with 
f7i } 2 the largest energy scale, all atoms from the most 
strongly interacting combination a = 1, 2 will form pairs, 
while those with a = 3 remain unpaired. Trion forma- 
tion is suppressed by the constraint. At low energies we 
obtain the effective model H « H e ^ = Hsf + Hbf + Hzb 
[12], with 

H 3B = -(2JV|^i, 2 |)^[26^ J -m 3 >3,i], (2) 

(id) 

where the lattice gradient is defined as Vi,\t 3 ^ = 
^3,i+e A — ^3,i , &\ is the unit vector in the primitive lat- 
tice directions, fh a ^ = &^&<r,i and H 3 b for the bound 
state is obtained by integrating out the fermions (1,2) 
in perturbation theory. The effective Hamiltonian H e fi 
provides us with an asymptotic strong-coupling picture 
in terms of a Bose- Fermi mixture, where fast fermions 
of the unpaired species with hopping J strongly interact 
with slow bosonic pairs, whose effective hopping param- 
eter is given by 4J 2 /|J7i, 2 | ~ J/10 <C J. The strong 
(off-site) boson-fermion repulsion J in Hbf is a direct 
consequence of the three-body constraint. This is the 
term that gives rise to phase separation - however, this 
scale is of the same order as the fermionic hopping and 
phase-separation is not generic for Bose-Fermi mixtures 
[13] , so the problem must be treated carefully. 
We assess the instability of the homogeneus system 



against phase separation by studying the compressibility 
matrix M aa > = d 2 E / dn a dn a > , where n ai n G i = n 3l m 3 
and E is the canonical energy for Eq. For low fill- 
ing and at mean-field level, we obtain E^ D \n 3 ^m 3 ) ~ 

* 2 mr2 {^) 2/D 4 /D+1 + Wy(4D\UiM + Wm 3 n 3 , 
where ilp = 27r,47r for D = 2,3. Consequently 
the system will be stable provided 1 < A = 
^{D/Slnf/CD-z^nl' 13 - 1 . For n 3 = 1/6, A = 

0.07(0.11) in 3D(2D), and thus the homogeneous system 
is unstable against phase separation. Note that there is a 
critical density below which the homogeneous mixture is 
stable in 3D (n w 10 -4 for this setup). In ID the system 
is not expected to phase separate on typical experimental 
scales, where the constraint helps instead to stabilize the 
homogeneous ACS phase [TO] . 

Since in strong-coupling T c = 0{ J 2 /\U la \) < J [14] . 
we expect phase separation (with a typical energy scale of 
the boson-fermion repulsion J) to be present well above 
the critical temperature for the disappearance of the SF 
phase. Finally, as the fermion-boson interaction is re- 
pulsive, the instability will result in spatially separated 
domains of fermions and bosons in order to minimize the 
bulk energy, as shown numerically below. 

Numerical approaches - In order to treat the system 
quantitatively in a wider range of parameters and go 
beyond the effective description above, we employ two 
complementary numerical approaches. Details on both 
methods are provided in [TO] . 

In 3D, we employ DMFT as a non-perturbative ap- 
proach where the (quasi-) momentum dependence of the 
self-energy is neglected, i.e. S(k, ui) = Z)(o;)[TO]. We 
generalized DMFT to the case with three species in or- 
der to describe ACS and trionic phase, and also inves- 
tigated the system at finite temperature. As a solver 
for DMFT, we used Exact Diagonalization [17]. Onsite 
thermodynamic observables, e.g. densities n a = (c^c a ), 
double occupancies d a ^> = (n a n a r) and the SF order pa- 
rameter P = (C1C2), are estimated directly as averages 
on the impurity site of the auxiliary Anderson model at 
self-consistency [TO] . 

The momentum distribution n a (k) and spectral prop- 
erties are obtained from the normal components of 
the single-particle Green function G a (k, iuj n ), which 
is a natural output in DMFT. Then, n a (k) = 
T^ n G a (k,iuj n )e~' lu;ri0 and the average kinetic energy 
per lattice site is given by K = jj- ^ k a £kn cr (k), where 
J\f is the number of lattice sites. The internal en- 
ergy per site E is then E = K + V pot: where V pot — 
J2 a <a' Ua,a'd a ,cT' is the mean potential energy per site. 

In 2D, we work with a VMC method, based on a 
strong-coupling expansion of Eq. [15 . We can de- 
termine the stability of the ACS phase with respect to 
trion formation and its instability towards phase sepa- 
ration by comparing trial wavefunctions for the ground 



state. For the normal phase we use 

\NFF) = JVzVdH J] O >> ( 3 ) 

where Vd is a fixed number projection, P 3 is a projector 
enforcing the three body constraint, and J is a Jastrow 
variational factor that can also include off-site trion cor- 
relations. For the homogeneous ACS phase we use 

\ACS) = JVsV D [] [« k + « k cl k>1 4 >2 ] J] 4', 3 |0>, 

(4) 

where and are the usual BCS functions. 

Since DMFT and VMC work respectively in the grand- 
canonical and canonical ensemble, they provide us with 
complementary descriptions of the system. This is par- 
ticularly relevant in the presence of phase separation. 

Numerical results at T = - We first treat the strong- 
coupling regime (A = 1), and characterize the ground 
state properties as a function of the overall density n for 
fixed interactions. In particular, we establish the insta- 
bility of the homogeneous ACS phase beyond the range of 
validity of the effective Hamiltonian approach and char- 
acterize quantitatively the equilibrium mixture. 

Within DMFT we find no values of the chemical po- 
tentials \i G such that the system is in a single homoge- 
neous phase with the required densities n a = n. This 
gives evidence for phase separation, in agreement with 
the general arguments above. This scenario, in which a 
mixture of phases has lower ground state energy E than 
a single homogeneous phase at T = 0, is signalled in the 
grandcanonical framework by the existence of multiple 
solutions with the same grand potential per lattice site 

= E - Y^a ^<y n <y for S iven I 1 *- 

We considered several possible scenarios to build the 
equilibrium mixture, combining phases with pairing in 
different channels [18] and also explicitly suppressing su- 
perfluidity. We found that the equilibrium mixture at 
T = consists of a fraction a of a two-species SF in- 
volving the species with the strongest interaction, 1-2, 



at densities 



SF 



n: 



SF 



n n , n 



SF 



0, and a fraction 



1 — a of a normal phase which accommodates the remain- 
ing unpaired species, i.e. = ni? = 0,n^ = n u . The 
energy of the mixture per lattice site is given by E mix = 
aE SF (n p ) + (1 — a)E N (n u ) and a = n/n p = 1 — n/n u . 
Note that a similar scenario (though with a different a) 
also applies for much smaller A, where the simple ana- 
lytical picture above does not apply. Moreover, we find 
that the densities in the two phases in equilibrium are 
very different, being n p ~ 0.75 and n u ~ 0.2 for A = 1 
and n = 0.16, since a (strongly- interacting) two-species 
SF is in equilibrium with a spin-polarized (and therefore 
non-interacting) Fermi gas. 

This scenario is confirmed via VMC in 2D, where we 
use an explicit Maxwell construction. We find that the 
energy E mix of the SF plus normal mixture described 
above is lower than the energy Eh om of the homoge- 
neous ACS solution for the same overall fixed densities 
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FIG. 2: (Color online) (a-c) VMC and DMFT results at T = 
0, for A = 1 as a function of n: a) Energy gain AE from 
phase-separation within VMC. b) Ground state energy E m ix 
of the equilibrium mixture, c) Fraction of superfluid phase a 
and superfluid order parameter P at equilibrium, d) Same as 
(c), but for n — 0.16 as a function of A. e) Same as (c), but 
for A = l,n = 0.16 at finite T. Energies are given in units 
with W = 1. 



n a = n. We plot AE = E hom - E mix > in Fig. [2^, 
where it is evident that the energy difference between 
the homogeneus solution and the equilibrium mixture is 
a very small fraction of E m i x shown in Fig. [2}d. The 
relative energy gain from phase separation is therefore 
very small, probably because the kinetic energy of the 
unpaired fermions and the boson-fermion repulsion are 
of the same order ~ J (see Eq. |2|. Once rescaled with 
the bandwidth W = 4DJ, the DMFT and VMC values 
for the energy E^f^ FT and E^f^ are in good agree- 
ment. The approximately linear scaling of E mix with the 
density (E mix « nU\^) occurs because the interaction 
energy in the superfluid is much larger than the kinetic 
energy in the superfluid or normal fluid, and ~ n P 
in the superfluid phase at strong-coupling. 

To further characterize the equilibrium mixture, 
we also evaluated the SF order parameter P = 

lim R ^ 00 P(R) = yj Ei(4,2 c !,i c i+R,i c i+R,2), where 
Nsf is the number of lattice sites in the SF domain. 
As shown in Fig. a increases monotonically with the 
density in both 2D and 3D. Conversely, P is maximum 
whenever n p = 0.5, similarly to a two-component SF. 
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FIG. 3: (Color online) (a,b) ^-Integrated momentum distri- 
bution n a (k x ,k y ) from DMFT for (a) a = 3 in the normal 
phase and (b) a = 1 (identical to a = 2) in the SF phase, (c) 
Raman response i2(ct>) for the normal (dotted blue line) and 
SF component (solid red line) of the mixture. All results are 
for A = 1, n = 0.16 and T = 0. 



In Fig. [2]i we show that the phase-separated scenario 
persists in the intermediate-coupling regime where the 
effective Hamiltonian description fails, rescaling the in- 
teractions with a factor A down to A = 0.3. 

Numerical results at finite temperature - To connect 
with experiments, we investigated the effect of finite 
temperature T within DMFT. For finite T the grand- 
canonical potential is given by Q = E — TS — ^ a ii G n G , 
where the entropy S is computed via the Maxwell re- 
lation drier /dT = dS/dficr. For increasing T the order 
parameter P for the superfluid component decreases and 
eventually vanishes at T = T c , as shown in Fig. [2k 



However, phase separation persists for T > T c , where 
both components are normal, as clearly seen from the 
mixing fraction <x As explained above, this effect is due 
to the loss- induced constraint and scale separation. This 
is in contrast to (i) strong coupling in the absence of a 
constraint, where the trionic phase is favored and (ii) the 
case of symmetric interactions U ajCr > = U, where phase 
separation is favoured due to a condensation energy gain 
in the ACS phase and disappears above the T c [T5] . 

Experimental signatures - Clear experimental signa- 
tures of this scenario appear in the quasi-momentum 
distribution of the mixture, accessible in state-selective 
time-of-flight (TOF) measurements. Within DMFT 
we estimated the z-integrated momentum distribution 
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(k) of the two phases in 3D at 



T = 0. For the normal phase of the unparied species 
(Fig. [3k), this coincides with the Fermi distribution at 
T = \Z = 1, where Z is the quasiparticle weight). De- 
spite smoothing due to the integration along the z axis, 
ft vanishes identically at large k. If the system were in 
a balanced ACS phase, the jump at the Fermi surface 
of the unpaired species would be strongly renormalized 
at strong coupling (Z < 1) and some particles would be 
promoted to momenta outside the Fermi surface. The SF 
component (Fig. [3]o) shows the weak dependence of 
on k as expected for a two-color SF in the BEC regime. 

Complementary signatures can be obtained using Ra- 
man spectroscopy [19 , which allows us to distinguish dif- 
ferent scenarios via the spectral gap. Within DMFT, 
the Raman response at zero transferred momentum 
R a (u)), where Co = uo + 11 — e Q out and s® ut is the en- 
ergy of the final state, is given at T = by R a (^) — 
^Ek k< "Ax(k,£ k - Q\ where C is a constant and 
A a (k,uj) = -l/7rImG a (k,cj + zO + ). As shown in Fig.J^, 
R(uj) for the SF component of the mixture shows a 
gap A « 3.2W = 0{U\2), while for the normal (non- 
interacting) component R(uj) oc 5(uj— ^3). For a balanced 
ACS, the Raman peak at uo = /i3 would be broadened due 
to the interactions and in a trionic phase every compo- 
nent would be gapped with A a = 0(U a , a +i + ^,^+2)- 
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